Turbulent Diamagnetism in Flowing Liquid Sodium 
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The nature of Ohm's law is examined in a turbulent flow of liquid sodium. A magnetic field is 
applied to the flowing sodium, and the resulting magnetic field is measured. The mean velocity 
field of the sodium is also measured in an identical-scale water model of the experiment. These two 
fields are used to determine the terms in Ohm's law, indicating the presence of currents driven by 
a turbulent electromotive force. These currents result in a diamagnetic effect, generating magnetic 
field in opposition to the dominant fields of the experiment. The magnitude of the fluctuation-driven 
magnetic field is comparable to that of the field induced by the sodium's mean flow. 
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The magnetic fields of stars and planets are generated 
by the motion of an electrically conducting plasma or 
liquid metal in the star or planet's interior The mo- 
tion of the fluid across an existing seed magnetic field 
generates a motional electromotive force (EMF) that 
drives currents, leading to a self-generated magnetic field. 
While large scale fiow may induce much of the mag- 
netic field, there is another means by which magnetic 
field may be generated: the interaction of velocity and 
magnetic field fiuctuations [5]. If the magnetic and ve- 
locity fields are separated into their mean (taken to be 
the temporally averaged value) and fluctuating compo- 
nents, B = (B) -I- B and V = (V) -t- V, then the mean 
electric current obeys an Ohm's law of the form 



J) = a ((E) + (V) X (B) + (V X B) 



(1) 



where (J) is the average current density, a the conduc- 
tivity of the fluid, and (E) the average electric field. The 
electric field plays a passive role in the scenario stud- 
ied here, merely maintaining V • (J) = 0, and as such 
will not be considered further. There are two significant 
source terms in Eq. 1: (V) x (B) represents the EMF 
associated with the mean part of the velocity and mag- 
netic fields, while (V x B) represents the EMF generated 
by the fluctuating part of the fields; velocity and mag- 
netic field fluctuations can interact coherently to generate 
mean currents. 

The possibility of such a turbulent EMF has long been 
recognized 0, 0, Much of its study has focused on 
those currents generated by helical velocity field fluctu- 
ations 0, 0|, and several non-simply-connected liquid- 
metal experiments have been constructed to mimic such 
helical flows 0,@|. However, gradients in the intensity 
of the fluctuations can also generate currents. These 
currents effectively expel magnetic flux from regions of 
high turbulence to low, resul ting in a diamagnetic ef- 
fect 0, i, [i3Jll H [li El Ell • Such flux expulsion 
may explain [16| the weak magnetic fleld at the center 
of the galactic core [17| relative to the core's external 



flux tubes [18|, as well as concentrations of large scale 
toroidal magnetic field at the base of the stellar convec- 
tion zone [ii,[23,[2H. However, since fluctuation-driven 
fields cannot yet be isolated in astrophysical or geophys- 
ical settings, whether turbulent EMFs play a significant 
role in the production of such magnetic fields remains an 
open question. 

In this Letter we present the spatial structure of such 
fiuctuation-driven currents, as measured in the Madison 
Dynamo Experiment. The net result is a strong turbu- 
lent diamagnetism, reducing the magnitude of the mag- 
netic field throughout the experiment. A previous paper 
showed that such currents must exist |22l |. but gave no 
information about their structure. Here, direct measure- 
ment of the experiment's magnetic and velocity fields is 
used to explicitly examine the structure of Ohm's law. 

The Madison Dynamo Experiment is a one-meter- 
diameter sphere of flowing liquid sodium [23| . Sodium 
is chosen as the working fluid for its high electrical con- 
ductivity. An axisymmetric mean velocity field is gener- 
ated within the sodium by a pair of counter-rotating im- 
pellers attached to shafts that enter the sphere through 
each pole (with the shafts defining the axis of symme- 
try). The flowing sodium is very turbulent, with a kinetic 
Reynolds number Re ^ 10^. The unconstrained geome- 
try of the experiment allows fluctuations in the velocity 
field to develop up to the scale of the device. The fluctu- 
ating velocity field generates a fluctuating magnetic field 
by advecting the mean magnetic field of the experiment. 

The nature of Ohm's law is explored by measuring the 
magnetic-field-dependent terms in Eq. [T] to determine if 
the induced magnetic field is due solely to the action 
of the mean flow. An approximately-uniform magnetic 
field is applied to the fiowing sodium using a pair of ex- 
ternal magnetic field coils, and the total magnetic field, 
applied plus induced, is measured (the applied field is suf- 
ficiently weak that the velocity field is unaffected). The 
collection of magnetic field data has been described pre- 
viously [2^ . Since both the mean velocity and applied 
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FIG. 1: Induced magnetic fields measured in the Madison Dynamo Experiment, (a) Mean measured poloidal and toroidal 
field values, scaled to the magnitude of the applied field (50 Gauss), for an impeller rotation rate of 1000 RPM. The fit (solid 
lines) represents values predicted by a spherical harmonic expansion fit to the data (diamonds). Error bars represent the RMS 
fiuctuation levels of the signals, (b) The reconstructed field, with the axis of symmetry oriented horizontally. Streamlines of 
the poloidal field are in the upper hemisphere and the toroidal field strength is in the lower hemisphere. Measurement positions 
are indicated with dots, and the impeller positions are indicated with rectangles. 



magnetic fields are axisymmetric, the mean induced field 
is also axisymmetric (all data presented in this Letter 
are axisymmetric). A spherical harmonic expansion of 
the induced internal magnetic field, in the poloidal and 
toroidal directions, is fitted to the most probable val- 
ues of the measured magnetic field. Since the sphere is 
simply connected, a toroidal magnetic field cannot be 
applied from outside the sphere; all measured toroidal 
magnetic field is due to electrical currents fiowing within 
the sodium. 

For an impeller rotation rate of 1000 revolutions per 
minute (RPM), the mean measured induced magnetic 
field ranges from 1.2 times the magnitude of the ap- 
plied field in the poloidal direction, to 1.1 times in the 
toroidal direction [Fig. [T]^a)]. The field reconstructed 
from the fit indicates that the external poloidal magnetic 
field is dominated by a dipolar component that opposes 
the externally-applied magnetic field [Fig.[ljb)]. The re- 
construction also demonstrates that the toroidal velocity 
field is effective at generating toroidal magnetic field from 
the applied poloidal magnetic field. 

To distinguish between (B)^-^^^^^^, the magnetic 
field induced by the mean flow (V x (B)^.^^^^^^ = 
Mo{J)(v>x(B) ^ /^ocr(V) X (B)), and (B)^vxb)' t^e 
magnetic field due to fluctuations (V x (B)^.^^3^ ^ 

/Ltocr(V x B)), the mean velocity field of the sodium must 
be known. To this end an identical-scale water model of 
the sodium experiment has been constructed . At the 
correct temperatures water and sodium have the same 



kinematic viscosity and similar densities. As a result 
the two fluids are nearly hydrodynamically identical and 
water can be used to model the flowing sodium. Like 
the sodium apparatus, the water model is a one-meter- 
diameter sphere in which impellers generate an axisym- 
metric mean velocity field. Stainless steel tubes, identical 
to those containing the internal Hall probes in the sodium 
experiment [l^], enter the fiow at the same seven loca- 
tions as in the sodium experiment. Unlike the sodium 
apparatus, the water model is outfitted with five windows 
that allow a Laser Doppler Velocimetry (LDV) system to 
directly measure the poloidal and toroidal components of 
the water's velocity field. A spherical harmonic expan- 
sion of the velocity field is fitted to the mean values of 
these measurements to determine the mean velocity field 
of the flowing sodium. For a impeller rotation rate of 
1000 RPM, the measured velocity field data yield a max- 
imum poloidal speed of 2.5 m/s and a maximum toroidal 
speed of 3.9 m/s [Fig. [DJa)]. As expected, the fiow is 
counter-rotating in the toroidal direction and the poloidal 
flow rolls inward at the equator and outward at the poles 
[Fig.Hb)]. 

Once the mean velocity field is known the magnetic 
field due to the mean flow interacting with the mean 
magnetic field [Fig. [Slja)] and the magnetic field driven 
by the fiuctuations [Fig.[3jb)] are calculated. Like the fit 
to the measured magnetic and velocity fields, this calcu- 
lation is done in a spherical harmonic basis, and involves 
calculating the electrical potential of the experiment as- 
suming V • (J) = 0. If no fiuctuation-driven currents 
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FIG. 2: Velocity fields measured by Laser Doppler Velocimetry in the water model of the Madison Dynamo Experiment, (a) 
Mean measured velocity field as a function of radial position, for an impeller rotation rate of 1000 RPM. (b) The reconstructed 
field. The poloidal fiow rolls inward at the equator and outward at the poles. Two toroidal cells rotate in opposing directions 
in each hemisphere. In both panels the conventions are as in Fig. [T] 



were present (i.e. the magnetic field is solely explained 
by (J){v)x(B) ~ (^) ^ (■^))' then the measured mean 
magnetic field [Fig. [TJ^b)] would be the same as the field 
due to the mean velocity field [Fig. [3)Ja)]. Inspection 
of these two figures reveals that this is not the case. 
The prominent external dipole component of the mea- 
sured poloidal magnetic field is completely absent from 
the magnetic field due to the mean velocity field, since the 
mean axisymmetric velocity field is incapable of produc- 
ing it [2^. Also, the magnitude of the measured toroidal 
magnetic field is significantly weaker than the toroidal 
field due to the mean velocity field interacting with the 
measured magnetic field. Clearly, the inductive action of 
the mean flow alone is insufficient to explain the mea- 
sured fields: fluctuations must be generating significant 
magnetic field. 

The magnetic field generated by the fluctuations has 
several prominent features [Fig. [3]^b)]. First, there is a 
dipole component that dominates the magnetic field out- 
side the sphere, in the direction opposite the applied field. 
This is the source of the measured external dipole mo- 
ment. Second, the toroidal field induced by the fluctua- 
tions is in opposition to the fleld induced by the mean ve- 
locity field interacting with the measured magnetic field. 
Thus, the field induced by the fluctuations is diamag- 
netic with respect to the dominant poloidal and toroidal 
fields within the experiment. The effect is important 
to the overall magnetic field of the experiment, as the 
fluctuation-driven field is a significant fraction of both the 
field induced by the mean flow, and the applied field. The 
strength of the toroidal diamagnetic field near the poles 



is about 50% of the magnitude of the field induced by 
the mean flow, while the fluctuation-driven dipole com- 
ponent is 20% of the magnitude of the applied field. The 
fluctuation-driven poloidal field reduces the total poloidal 
flux through the equatorial plane by 10%. 

Characterization of the fluctuations which lead to this 
turbulent EMF requires direct measurement of the fluc- 
tuating components of V and B. Since such measure- 
ments are not currently available, it is not yet possible 
to determine whether the observed diamagnetism is due 
to gradients in fluctuation levels, helical turbulence, or 
some other effect. Since the fluctuation levels of the ex- 
periment are greatest near the impellers, and weakest 
near the sphere's surface, flux expulsion due to gradients 
in turbulence levels is a natural candidate to explain the 
observed diamagnetism, though it is not the only possible 
explanation. 

We note that the structure of the diamagnetic fleld is 
qualitatively similar to the fluctuation-driven fleld pre- 
dicted by recent numerical simulations of the experi- 
ment '25']. These simulations predict both the presence 
of the dipole moment and the overall weakening of the 
measured toroidal fleld. However, the magnitude of the 
fluctuation-driven fleld presented here is a factor of five 
larger than that predicted by the simulations. The rea- 
son for this discrepancy is not currently known, though 
it may be related to differences in fluid forcing or the 
simulation's large value of magnetic Prandtl number. 

In summary, mean magnetic and velocity field mea- 
surements have been used to determine the structure 
of fluctuation-driven currents in the Madison Dynamo 
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FIG. 3: Magnetic field due to the terms in Ohm's law. (a) The magnetic field, (B) (y)x{'B) calculated from the mean measured 
velocity field interacting with the measured mean magnetic field (the sum of the induced and the applied fields). Note the lack 
of induced external dipole moment, (b) The mean magnetic field, due to the EMF associated with the fiuctuations. 

The poloidal external field displays the measured dipole component, in opposition to the applied field. The toroidal field is in 
opposition to the field induced by the mean fiow of sodium. 



Experiment. These currents lead to a magnetic field in 
opposition to the experiment's dominant magnetic field. 
This is the first observation of turbulent diamagnetism 
in a laboratory setting. Analysis of the nature of the 
fluctuations which leads to this turbulent EMF is ongo- 
ing, and future work will focus on further understanding 
the localized regions of current generation. Given that 
the observed effect is diamagnetic, it indicates that such 
fields could be a hindrance to the development of mag- 
netically self-exciting geophysical and astrophysical sys- 
tems, as well as simply-connected liquid-metal dynamo 
experiments. 
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